STRESS  ANALYSIS,    BUCKLING  ANALYSIS   AND 
OPTIMUM   PROPORTIONS   OF  AN   ISOTROPIC 
WEB-STIFFENED  SANDWICH   CYLINDRICAL 
SHELL  UNDER  HYDROSTATIC    PRESSURE 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Department  of  Naval  Architecture  and  Marine  Engineering 


STRESS  ANALYSIS,  BUCKLING  ANALYSIS  AND  OPTIMUM  PROPORTIONS 
OF  AN  ISOTROPIC  WEB-STIFFENED  SANDWICH  CYLINDRICAL  SHELL 

UNDER  HYDROSTATIC  PRESSURE 

by 

FRANK  ROBERT  KRONER,  JR. 

B.S.,  United  States  Naval  Academy 
(1961) 


SUBMITTED  IN  PARTIAL  FULFILLMENT 
OF  THE  REQUIREMENTS  FOR  THE 
DEGREES  OF  NAVAL  ENGINEER 
AND  MASTER  OF  SCIENCE  IN 
OCEAN  ENGINEERING 

at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
May,  19  7  0 


>     * 


4L  P0STGRA3&ATS  8C 

koiiiTEBEY,,  CALIF.. 93.940^" 

2, 

STRESS  ANALYSIS,  BUCKLING  ANALYSIS  AND  OPTIMUM  PROPORTIONS 
OF  AN  ISOTROPIC  WEB-STIFFENED-  SANDWICH  CYLINDRICAL  SHELL 

UNDER  HYDROSTATIC  PRESSURE 
by 
Frank  Robert  Kroner,  Jr. 

Submitted  to  the  Department  of  Naval  Architecture  and  Marine 
Engineering  on  May  21,  1970,  in  partial  fulfillment  of  the 
requirements  for  the  degree  of  Naval  Engineer  in  Naval 
Architecture  and  Marine  Engineering  and  Master  of  Science  in 
Ocean  Engineering. 

ABSTRACT 


Web-stiffened  sandwich  cylindrical  shells  consisting  of 
two  concentric  cylindrical  shells  joined  and  stiffened  by 
annular  rings  (webs)  were  considered  in  the  investigation. 
The  purpose  of  this  investigation  was  to  study  the  stress 
distribution  and  buckling  modes,  and  to  determine  the  scant- 
lings and  configuration  of  optimum  proportions. 

The  scantlings  and  configuration  of  optimum  proportions 
was  defined  as  that  with  the  lowest  weight  to  displacement 
ratio  at  a  given  loading  of  hydrostatic  pressure. 

The  scantlings  (outer  shell  thickness,  inner  shell 
thickness,  web  thickness,  web  depth  (separation  of  shells) 
and  configuration  (web  spacing)  were  systematically  and 
independently  varied  and  the  effect  upon  failure  mode  pres- 
sures, and  the  weight  to"  displacement  ratio  v/as  determined. 
The  results  were  plotted  graphically. 

The  form  of  stress  and  buckling  analysis  results  follow 
closely  those  of  a  ring-stiffened  cylinder.   Geometric  sim- 
ilitude was  determined  to  exist.   The  scantlings  and  configu- 
ration of  optimum  proportions  were  found  to  be  a  function  of 
the  loadings. 

Web-core  sandwich  configurations  possess  structural 
efficiencies  on  the  order  of  20%  higher  than  those  of  con- 
ventional ring-stiffened  construction  for  loadings  on  the 
order  of  1000  psi  (2225  feet)  and  these  diminish  with  depth 
to  about  5%  higher  for  loadings  on  the  order  of  4  000  psi 
(8900  feet) . 

The  real  advantage  of  the  web-stiffened  sandwich  con- 
struction lies  in  the  use  of  thinner  plating  with  all  its 
superior  metallurgical  and  fabrication  qualities. 

Thesis  Supervisor:   Professor  J.  H.  Evans 

Title:  Professor  of  Naval  Architecture 
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b  Web  thickness 

c  (W  '  1 

C  Web  constant 


D  Bending,  rigidity,  Et  /12(l-v  ) 

D  Mean  diameter  of  shells 

M 

d  RRO-RRI 

E  Young ' s  modulus 

E  Secant  modulus 

s 

E  Tangent  modulus 

f  Stress  ratio,  a  /a, 

g  ,  g.  Edge  coefficients  for  outside  and 

inside  shells 

a  ,    a    . ,    a .     ,    a. .  Lame  deflection  coefficients  for  the 
^oo   ^01   ^io'  ^11 

(M,N,P,Q)         annular  webs 
h  R0-Ri 

k.  Stress  intensity  per  pound  of  pressure 

(Kx  =  ax/p) 

L  Unsupported  length  of  cylinder 

L  Length  between  bulkheads 

m,  n,  k  Integers 

aj  ,  iNj  Longitudinal  and  circumferential 

x'  $ 

forces  per  unit  length 

x,  ({>,  r  Longitudinal,  circumferential,  and 

radial  coordinates 
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B.  ,  F.  /  G.        Hyperbolic  and  trigonometric  stress 

parameters  for  outside  shell 

B  .  ,_ ,  F._.  G.T      Hyperbolic  and  trigonometric  stress 
il'   il    ll       J^  ^ 

parameters  for  inside  shell 

Longitudinal  moment  per  unit  length 

External  hydrostatic  pressure 

Force  coefficients  for  outside  shell 

Force  coefficients  for  inside  shell 

Total  radial  forces  on  the  outer  and 
inner  surface  of  webs  per  unit  length 

Longitudinal  forces  on  outside  and 
inside  shells  per  unit  length 

Radial  forces  on  outside  and  inside 
shells  per  unit  length 

Same  quantity  used  in  general  theory 
and  computer  programs 

m 

Radial  force  per  unit  length 
Radius 

Mean  radius  of  outside  and  inside 

shells 

Extreme  radius  of  outer  fiber  of  out- 
side  and  inside  shells 

RRI/RRO 

Mean  distance  between  webs 

Thickness  of  shell 

Outside  shell  thickness 

Inside  shell  thickness 
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u,  v,  w  Longitudinal  circumferential  and 

radial  deflections 

WQS  ,  tfj_s  Radial  deflection  of  outside  and 

inside  shell 

WQW,  Wj  Radial  deflection  of  the  outer  and 

inner  surfaces  of  the  web 

Z  Radial  bending  rigidity,  E/(l-v2) 

e  ,  e,/  e  Longitudinal,  circumferential,  and 

x'   <J>   r  r  ' 

radial  strains 

a  ,  aw  a  Longitudinal,  circumferential,  and 

radial  stresses 

Qq/6^  Shell  flexibility  parameters  for  out- 

side and  inside  shells, 
[3(l-v2)]  Vh     [L/(Rt)  l/2] 

ty  Bending  stress  coefficient, 

i(l-v2)/3]  V2 

v  Poisson's  ratio  for  isotropic 

materials 

a.  Stress  intensity,  Huber-Hencky , 


Von  Mises,  [  .'a?  +  oi    ~  °  \a  2*  ^  ^ 


a  Yield  stress 

y 

a.  . VT  (a.r,T)  Stress  location 

i  =  x,$,r      direction 

j  =  F,M        shell  (web-shell  juncture,  midbay) 
K  =  0,I,W      outside  shell,  inside  shell,  web 
L  =  0,I,M      fiber  (outer,  inner,  membrane) 

(jiKL)  i.e.,  MXOI  =  midbay,  longitudinal,  outside 

shell,  inner  fiber 

P  Yield  failure  at  shell  midbay 


PF  Yield  failure  at  web-shell  juncture 

PWS  Yield  failure  at  web  outer  fiber 

PPC1  Axisymmetric  shell  buckling 

PPC2  Asymmetric  shell  buckling 

PPC3  General  instability 

PPC4  Web,  ring  buckling 

PPC5  Web,  edge  compression  buckling 


INTRODUCTION 

One  of  the  more  promising  structural  concepts  for  the 
design  and  fabrication  of  cylindrical  pressure  hull  struc- 
tures is  the  sandwich  concept.   Structural  engineers  in  the 
aircraft  industry  nave  long  recognized  and  taken  advantage 
of  the  favorable  strength-weight  characteristics  of 
sandwich-type  construction.   In  studying  the  literature, 
however,  it  nas  been  found  that  the  loading  conditions 
encountered  in  tnese  applications  nave  dictated  sandwich 
structural  arrangements  that  would  be  of  no  direct  use  in 
the  design  of  pressure  hulls  for  submersibles . 

The  demands  of  hydrostatic  pressure  loading  are  such 
that  in  order  to  exploit  tne  sandwich  concept  for  pressure 
hull  construction  a  core  comprised  of  elements  which  are 
compression   resistant  as  well  as  shear  resistant  is  sought, 
The  major  difference  is  that  membrane  loads  are  predominant 
in  hydrospace  applications,  wnereas  in  aerospace,  the  bend- 
ing and  shear  type  loads  are  of  prime  concern. 

Experimental  programs  and  analytical  studies  have  been 
going  on  concurrently  in  order  to  develop  rational  formulas 
based  on  thin-shell  theory  for  predicting  the  static  struc- 
tural response  of  these  types  of  structures.   In  reference 
10  Pulos  presents  an  analysis  of  the  axisymmetric  elastic 


deformations  and  stresses  in  a  web-stiffened  sandwich 
cylinder  under  hydrostatic  pressure.   Raetz  presents  a 
similar  analysis  for  the  toroidal  tube-stiffened  sandwich 
cylinder.     In  references  12  and  14  Nott  presents  a  sim- 
plified stress  ana  strain  analysis. 

Nott ' s  equations  were  used  and  programmed  for  the 
stress  analysis  of  the  web-stiffened  structure. 

In  reference  22  the  axisymmetric  and  asymmetric  buck- 
ling equations  are  developed.   In  references  15  and  23  an 
expression  for  the  general  instability  is  determined. 

The  above  equations  are  utilized,  with  others  derived 
for  web  stability,  to  study  the  parametric  nature  of  web- 
core  sandwich  structure  and  to  determine  optimum  propor- 
tions . 

Optimum  design  must  be  based  on  rational  considera- 
tions of  inelastic  behavior,  however,  and  not  on  the  "one 
horse  shay"  concept  based  on  elastic  considerations  of  in- 
stability.  The  "ignorance  factors"  can  then  be  representa- 
tive of  the  variability  introduced  by  certain  intangibles 
which  are  not  easily  considered  in  a  theory. 

Some  of  the  intangibles  which  influence  static  strength 
and  complicate  the  problem  so  that  appropriate  design  for- 
mulas cannot  be  derived  on  purely  theoretical  grounds  are 
those  inherent  in  the  fabrication  process  itself: 
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(1)  Initial  stress. 

(2)  Residual  stress. 

(3)  Imperfect  circularity. 

(4)  Non-isotropic  or  non-homogeneity  of  material. 

(5)  Actual  boundary  conditions  at  stiffening  rings. 

Therefore,  in  view  of  these  complications,  the  devel- 
opment of  satisfactory  design  criteria  must  first  be  pre- 
dicted on  rigorous  matnematical  theory  with  its  concomitant 
idealizations,  and  then,  empirical  factors  derived  from 
test  data  can  be  introduced  to  "adjust"  the  theories  to 
take  account  of  these  many  variables  which  could  not  or 
were  not  considered. 

Thus  an  inelastic  analysis  of  all  instability  modes  is 
used. 
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GENERAL  THEORY 

1.  Stress  Analysis 

A  theoretical  analysis  of  the  ax i symmetric  elastic 
deformations  and  stresses  in  a  web-stiffened  sandwich 
cylindrical  shell  structure  under  external  hydrostatic 
pressure  is  presented  in  Appendix  A.   The  solution  is  based 
on  the  use  of  edge  coefficients  for  plate  and  shell  ele- 
ments of  finite  length,  and  includes  the  computation  of 
the  edge  forces  and  moments  arising  at  the  common  junc- 
tures of  these  elements. 

Equations  are  given  for  computing  numerically  the 
longitudinal  and  circumferential  stresses  in  the  two 
coaxial  cylindrical  shells  and  the  radial  and  tangential 
stresses  in  the  web  stiff eners  between  the  two  shells. 

To  evaluate  the  structural  strength  of  cylindrical 
sandwich  shells,  the  locations  and  magnitudes  of  maximum 
stresses  must  be  determined.   Rigorous  analyses  for  stress 
distributions  in  sandwich  shells  loaded  under  external 
hydrostatic  pressure  were  carried  out  by  Pulos    and 
Raetz   .   These  analyses  illustrate  that  the  maximum 
stresses  occur  in  the  inner  and  outer  shells  at  locations 
next  to  the  annular  webs  and  midway  between  the  two  webs 
that  separate  the  shells.   Therefore,  the  stresses  in  the 
inside  and  outside  shells  at  these  two  locations  are  of 
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interest  for  the  purpose  of  structural  design  and 
evaluation. 

The  procedure  and  nomenclature  for  programming  the 

12. 14 
stress  analysis  is  that  of  Nott. 

a.  NOMENCLATURE 
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b.  ASSUMPTIONS 

From  symmetry  considerations  it  is  seen  that  the 
edges  of  the  web  stiffener,  circular  annulus ,  do  not  under- 
go any  rotation.   This  stems  from  the  fact  that  a  horizontal 
tangent  or  zero-slope  condition  is  assumed  to  exist  at  the 
junctions  of  the  webs  with  the  two  cylindrical  shells. 
This  assumption  implies  that  the  edge  moments  on  each  shell 
at  the  shell-web  juncture  balance  each  other,  so  that  there 
are  no  net  moments  -co  be  resisted  by  the  web.   Further,  it 
is  assumed  that  the  web  elements  do  not  take  any  axial 
force  due  to  the  axial  pressure,  but  that  this  is  all  re- 
sisted by  the  cylindrical  shells. 

Also  inherent  in  the  analysis  (see  Appendix  A)  is  the 
neglect  of  the  beam-column  effect  due  to  the  axial  portion 
of  the  nydrostatic  pressure. 
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dw  +  lw  +   Etr=p_vN 

dx4  x   dx2         R2  R     x 


neglecting   the   beam-column   effect   becomes 


D|!w  +   EtWssp_vN 
dx*         R2  R      x 


Thus  the  analysis  of  the  web  stiffeners  reduces  to 
that  of  a  circular  annulus  subjected  to  axisymmetric  in- 
plane  radial  forces  on  both  its  inner  and  outer  boundaries. 

On  the  basis  of  these  assumptions,  it  is  only  necessary 
to  derive  edge  coefficients  for  an  annulus  undergoing  radial 
deflection. 

The  functions  A,    N,  P,  and  Q  are  the  Lame  deflection 
coefficients  for  the  annular  webs  and  can  be  expressed 

Rw     , ,     .       d    ,      , 2^w  d  _    . 

M  ■  b"  (1  +  2#    (—  +  2R^  -  2v) 

%    M  d    s      ,  2Rw  d  „. , 

N   "    -5-    U   +    2^»     (—  +    2^  -    2) 

Rw  d  2RW  d 

B^  d  2RW  d 

and  the  functions  g  ,  g .  are  the  edge  coefficients  from  the 
outside  and  inside  shells: 
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Thus,  from  Appendix  A,  the  axial  and  transverse  forces 

(V  ,  V.,  H   and  H . ) ,    which  act  at  the  intersections  of  the 
o    l    o      1  ' 

annular  webs  and  cylinders  can  be  determined. 
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And  from  these  forces  the  maximum  stresses  in  the 
shells  can  be  determined. 
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F20  =  9o 
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38_      sinh  6   -  sin  6 
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/3  (I-v2)   cosh  6   -  cos  6 
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a  =  /3(i-v2) 


Sr   t 
o  o 


and  likewise  for  B-^,  B^-,  .../  F,  .  ,  9  .  by  substituting  the 
inside  I  for  the  outside  0. 

d.  WEB  STRESSES 

From  Appendix  A  the  maximum  stress  in  the  web 

[ (-HHI (RR) 2+HH0)  -  (-HHI (RR) +HH0 (RR) ) (RR) ] 


a    =  -p 
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a^  =  ~P 


b/2  (l-(RR)2) 
[ (-HHI(RR) 2+HH0)  +  (-HHI (RR) +HH0 (RR) ) (RR) ] 


*WO  b/2  (l-(RR)2) 
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2 .  Buckling  Analysis  and  Failure  Modes 

The  axisymmetric  and  asymmetric  (lobar)  buckling 

;  4  8 

equations  developed  by  Lunchick   and  Reynolds   for  ring- 
stiffened  cylinders  are  modified  in  Appendix  3  for  applica- 
tion to  web-stiffened  sandwich  cylinders.   They  are  eval- 
uated for  both  the  inner  and  outer  shells.   Reference  22 
indicates  good  agreement  between  theoretical  and  experi- 
mental collapse  pressures. 

The  general  instability  equation  used  has  found  con- 

.  .    ,  .  •  .     ■*•«■*■    15,23 

sideraDle  experimental  verification. 

The  two  modes  of  web  buckling  are  first,  a  variation 
of  the  Foppl  or  Levy  ring  buckling  (effectively  TOKUGAWA) 
and  second,  the  buckling  cf  an  annulus  in  biaxial  edge 
compression  (see  Appendix  B) . 

a.  FAILURE  MODES 
1.  STRENGTH 

Using  the  Huber-Hencky-Von  Mises  criterion 


a  •    =      /a  i   +   a 2  -    o~~a~ 
i  $  x  cf)   x 


ai  "   p    ^Kj   +   k*    -   K^KX 


a 

p-  y 


K$   +   Kx    ~    KcpKx 


evaluated  at  the  web-shell  interface  (PF) ,  midbay  (P) ,  and 
outer  fiber  of  web  (PWS) . 
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2.  AXISYMMETRIC   BUCKLING 
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3.  ASYMMETRIC   BUCKLING 
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fE  E  VK* 


4.  GENERAL    INSTABILITY 


PPC3  = 


/E_ • E  •  • 
T   s 


*T 


n2-l+^-]  (n2+X2)  jl-3/4 


1- 


E  >| 

s 


X 


(n2+A2)2 


(n2-l)  Em  I 


T   EFF 


R3L 


where 


22 


A  = 


7TR 

LB 


ana 


"eff  for 


^  b  «- 


A 


t_  =  T   +  T. 
T     o     1 


R  =  R 


w 


5.  WEB    BUCKLING 


a. 4     Ring  Buckling 


PPC4  = 


24  E 

A*    £'T   eff 

dIl   (l-v2) 

m 


b.  Edge  Compression 


47T2Er 


PPC5    = 


12     (1-VZ)      (K^+4^)  { 


h 


23 


eff " 


k-L/2-*| 


T 

d 


1 

* 

-» 

€-b 

♦  , 

1 

bd3         LT^         LT3  /-d         T    12 

Xt  =   ~T2~  +   ~TT  +   TT  +   ToL      2    ~   T° 


T-  >2 


+    T.L     f|-^ 
1        I2  2  . 


A  =     (TL+T.L+    db) 
O  1 


y   =    (Y         -   d/2) 

2  eg  / 


eg 


ToL    (d   -  ^)    +   db    (|)    +   T.L(^) 

(T    L   +    db    +    T.L) 
o  i 


Xeff  =   ^  "  A* 
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AXISlTKttETRlC   SUCKLING-  (Bt'TWEEN   WE&S) 


ASYMNETRlC(LOBAR)BUCKLlNG(B&TWt£N  WEBS) 


GENERAL    INSTABILITY    ( CONCURRENT  FAU.U  RE 
OF   SHELLS  AND  WEBS) 


FAILURE  MODES  (SHELL) 
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W£b(  RlNG)SUCKLiN<i 


))  ))  » 


))  ))  )) 


WEB  (COMPRESSION)  0UCKLIWG 


FAILURE  IXODES  (WEB) 
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PROCEDURE 

The  buckling  pressure  is  a  function  of  the  cylinder 
geometry  and  the  secant  and  tangent  moduli  as  determined 
from  a  stress-strain  intensity  diagram  for  the  shell 
material . 

Before  the  buckling  equations  can  be  used,  E   and  E™ 

o  J. 

must  be  related  to  tne  applied  pressure.   The  secant  and 
tangent  moduli  are  defined 


s     i'  1 


ET  =  dc^/de. 


and  are  shown  graphically 
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Hence  E   and  Em  are  readilv  determined  from  the  stress' 

s       T  ■* 

strain  curve  (see  Appendix  C) .   Having  determined  E  /E  and 
ET/E  as  functions  of  a. ,    they  must  be  applied  to  the  hydro- 
staticaily  loaded  cylindrical  shell.   By  expressing  hydro- 
static pressure  in  terms  of  the  stress  intensity,  a  rela- 
tionship between  E  ,  E^,  and  pressure  will  be  established. 

O        J. 


<T\ 


(1)  Pc  =  f  (geometry,  Eg ,  ET) 
a  buckling  equation 

( 2 )  a .    =  /a?  +  at.    -   a   a 


1 


12 


p/K2  +  K2  - 


"1 


1"2 


a  stress  intensity  (Huber-Hencky-Von  Mises) 
equation 


The  intersection  of  these  two  curves  will  predict  fail- 
ure by  inelastic  buckling  if  below  the  yield  strength. 

Thus  the  general  solution  technique  is  to  (1)  solve 
for  the  state  of  stress,  and  consequently  a  relationship 
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between  J-  and  p,  and  then  (b)  solve  the  buckling  modes  as 
a  function  of  o. (E  ,ET) .   The  intersection  of  the  curves 
being  the  critical  pressure  (see  lb) . 


1.  Solution  Technique. 
a.  STRESS   ANALYSIS 
INPUT  DATA: 

Ro'  RI'  To'  TI'  b'  L'  V/  p 
Let  p  =  1.0  psi 

(1 

(2 

Formulate  g  ,    g. 
3o    1 


(3 
(4 
(5 
(6 
(7 
(8 
(9 
(10 
(11 


Formulate  3  ,  6- 

Formulate  F1Q,  F.^,  F2Q,  F21  ... 


Formulate  R  ,  d 
w 

Formulate  M,  N,  P,  Q 

Formulate  V  =>  V. 
o    i 

Solve  simultaneous  equations  for  H  , H . 
Formulate  HHO,  HHI,  RR 

Formulate  a   ,  a  . 

uo   ui 

Formulate  B1Q/  Bin,  B2o'  S21  '** 
Solve  for  the  states  of  stress 


b.  FAILURE   MODES 

(1)  INELASTIC  MATERIALS 
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PPC3 

w             ?y 

PPC1 

\y^  i 

PPC2 

V           J*?      i 

JR     \  ' 

jo              \       V 

\  V  \ 

\   \    \ 

!\ 

i 

<K 


(2)  ELASTIC  -  PERFECTLY  PLASTIC  MATERIAL 


ll 
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c.  OPTIMUM   PROPORTIONS 

A  simple  optimization  program  was  used  to  deter- 
mine optimum  proportions.   As  indicated  in  the  flow 
diagram,  a  comparison  is  made  between  the  design  col- 
lapse pressure  (PCO)  and  all  the  failure  mode  pressures 
(P.).   The  failure  mode  pressures  consist  of: 

P  =  yield  at  midbay 
PF  5  yield  at  web-shell  junction 
PWS  =  yield  in  web 
PPC1  =  axisymmetric  buckling,  each  shell 
PPC2  =   asymmetric  buckling,  each  shell 
PPC3  =  general  instability 
PPC4  =  web  instability  (ring) 
PPC5  =  web  buckling  (edge  compression) 


NPUTS 
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PROGRAM 


V 


GEOMETRY 
GENERATOR 


(no-go) 


COM  PRR^SOM 

(    fCOLtRPSc  'L   J 


(&o) 


OBJECT  IV  E 
PUN  C\OU 

(vieiGHT/^uoYRMCY) 


STRESS 

ANALYSIS 

SUBROUTINE 


FAILURE 
IXODES 
SUBROUTINE 
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RESULTS 


1.  Geometric  similitude  was  determined  to  exist.   The 
stress  analysis  and  failure  mode  programs  output  are  in 
Appendix  D  for  three  geometrically  similar  structures. 
Graph  1  shows  the  failure  modes  for  these  three  geometric- 
ally similar  structures  made  from  HY-160  steel.   The  effect 
of  the  inelastic  range  of  this  material  on  the  buckling 
modes  is  indicated.   The  example  has  underdeveloped  webs 

(PWS,  PPC4) . 

2.  Maximum  stress  states  were  found  to  be  located  as 
shown  below. 


wo 
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5 


3.    The  results  are  plotted  in  graphical  form.   The  results 
are  divided  basically  to  show  (a)  the  failure  modes  of  the 
shells,  (b)  the  failure  modes  of  the  webs,  and  (c)  optimum 
proportions  and  weight/displacement  results.   The  presenta- 
tion of  the  results  in  this  form  was  (1)  to  facilitate 
comparison  with  ring-stiffened  cylinders,  and  (2)  to 
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show  the  relationship  of  the  various  failure  modes. 
The  form  of  the  results  are  very  similar  to  those  of  ring- 
stiffened  cylinders. 

4.  Graphs  2  through  7  are  sample  results  of  the  stress 
and  buckling  analysis  for  the  shell  portion  of  the  web- 
stiffened  sandwich  cylinder  using  HY-30  steel. 

5.  Graphs  8  through  11  are  sample  results  for  the  web 
portion  non-dimensionalized  by  PXT_ 

where 

m    2(TQ  +  Tj)  gy 
PND  ^ 

The  material  was  HY-80  steel. 

6.  Graphs  12  and  13  show  the  effect  of  yield  strength  and 
the  modulus  of  elasticity  on  shell  failure  modes  respect- 
ively. 

7.  Graphs  14  and  15  and  Tables  1  and  2  are  concerned  witn 
the  optimization  results.   The  optimum  proportions  for  ring- 
stiffened  cylinders  shown  in  Graph  15  were  based  upon  using 
92,  92A  criteria  for  shell  yield  (outer  fibers  reaching 
yield  stress) . 

8.  For  a  given  hydrostatic  loading  and  hull  diameter  the 
shell  thicknesses  for  the  web-stiffened  sandwich  are  approx- 
imately half  that  of  a  ring-stiffened  cylinder. 


GRAPH  1 
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Table    1 
OPTIxMUM    SCANTLINGS 


HY-30 


P   , ,   h/R   h/s    h/t 
coll    '  w    '  '     s 


b/t 
s 

t  /D 
o' 

.0019 

t./D 
l 

L/D 

W/D 

.538 

.0015 

.020 

.137* 

.696 

.003 

.00275 

.030 

.239 

.600 

.005 

.005 

.025 

.436 

500  .058  1.43  17.65 

1000  .095  1.56  17.39 

2000  .113  2.14  12.00 

5000  .148  2.16    6.15   .538   .014    .014     .030   1.061 


HY-150 

P   -.   h/R    h/s    h/t     b/t    t  /D     t./D    L/D     W/D 
coll     w  s      s    o       1 


500  .077  2.50  39.02  .976  .0011  .0010  .015  .125 

1000  .077  1.87  26.67  .933  .0016  .0014  .020  .155 

2000  .077  1.73  14.29  .893  .0028  .0028  .020  .258 

5000  .113  2.14  9.23  .462  .008  .005  .025  .522 


This  value  seemed  high  and  thus  a  scatter  check  was  con- 
ducted about  this  point.   Results  of  this  check  are  in 
Table  2. 
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Table  2 
OPTIMUM  SCATTER  ABOUT 


P   nn  =  500  psi 
coll 


Pcoll   M°de   h/RN   fco/b    ti/D    L/D    b/ts   W/D    W/D1 


508  PF  .058  .0014  .0014  .025  .786  .124  .122 
o 

554  PWS  .058  ■.0014  .0014  .020  .643  .128  .115 

524  PWS  .058  .0015  .0014  .025  .759  .127  .121 

581  PWS  .058  .0016  .0014  .020  .600  .134  .115 

647  PPC4  .058  .0018  .0014  .020  .625  .145  .112 

677  PPC4  .058  .0018  .0015  .020  .667  .153  .113 


*                     500 
Normalized  W/D  by  (p -)  . 

coll 


66 

DISCUSSION  OF  RESULTS 

1.  High  strength  materials  having  inelastic  zones  are 
more  prone  to  failure  by  buckling  than  elastic  materials 
because  of  the  reduction  of  the  value  of  E  above  the  pro- 
portional limit. 

2.  Graphs  2  through  5  indicate  that  close  spacing  (small 
L)  of  the  annular  rings  (webs)  is  required  in  order  to 
utilize  the  benefits  of  this  type  of  structure.   This  fol- 
lows from  the  more  uniform  stressing  in  this  configuration. 

3.  Maximum  stresses  in  the  shell  occur  in  the  outer  shell 
(T   =  T.)  when  web  spacing  (L)  is  large.   However,  when  the 
spacing  is  small  the  stress  in  the  inner  shell  can  exceed 
the  outer.   This  occurs  when  the  web  thickness  and  web 
depth  is  insufficient  to  carry  the  load  (see  Graphs  6  and 
7) .   Graph  6  suggests  that  a  more  evenly  distributed  stress 
can  be  obtained  by  making  T   >  T.  for  thin  shells.   Optimi- 
zation bears  this  out  (Tab+e  1) . 

4.  Graphs  8  through  11  indicate  that  web  strength  and 
buckling  determine  web  depth  and  thickness  requirements  and 
have  a  critical  effect  on  the  weight/displacement  ratio 

(PWS  =  f(bh),  PPC4  =  f(bh3),  PPC5  =  f(b/h)2).- 
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5.  Shell  buckling  modes  (PPC1,  PPC2)  become  prevalent 
only  in  (a)  thin  shells  --  either  shallow  depth  load  design 
or  high  strength  materials,  or  (b)  materials  with  low 
moduli  of  elasticity  (E)  such  as  titanium  or  aluminum,  or 

(c)  in  the  ineiasric  range  of  high  strength  materials. 

6.  Axisymmetric  failure  mode  (PPCi)  is  prevalent  only  in 
very  thin  shells  (t/D   <  .0015). 

7.  Higher  strength  materials  and  materials  with  low 
moduli  of  elasticity  should  have  smaller  web  spacing  in  the 
optimum  proportions  (PPC2  shifts  to  the  left  --  thinner 
shells  for  the  high  strength  materials;  smaller  E  for  those 
with  low  modulus  of  elasticity) . 

8.  Table  1  and  Graph  14  indicate  that  web  depth  increases 
with  depth  for  the  optimum  scantlings.   This  is  seen  from 
the  need  to  increase  the  web  size  with  loading.   Also 

h/t  decreases  with  depth  indicating  t  is  increasing 
faster  than  web  depth,  h.  The  parameters  h/s  and  b/t 
do  not  change  appreciably. 

9.  Optimization  scatter  about  -che  optimum  scantlings 
indicates  more  than  one  optimum  configuration  (Table  2) . 

10.  Graph  15  shows  that  there  is  an  advantage  of  this 
structure  over  the  ring-stiffened  cylinder  of  some  order. 
Further  examination  is  required  using  membrane  yield  in 
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place  of  yielding  at  the  outer  fibers  to  compare  to  the 
optimized  results  of  ring-stiffened  cylinders  using 
Lunchick's  formulation  for  shell  yield.   Also,  comparison 
with  experimental  failures  is  required  to  verify  the 
formulation  used  in  this  investigation  or  to  justify  the 
use  of  membrane  stresses. 

11.  The  advantages  of  thinner  shells  (homogeneity  of 
material,  greater  ductility,  better  notch  toughness,  and 
easier  to  form  or  weld)  compared  to  those  of  the  ring- 
stiffened  cylinder  make  the  web-core  sandwich  cylindrical 
hull  very  appealing. 

12.  The  results  set  forth  are  thought  to  be  conservative. 
Yield  (P,  PF ,  PWS)  at  the  outer  fibers  was  the  criteria  of 
yield  failure.   If  membrane  stresses  were  permitted  to 
reach  yielding  more  dramatic  results  would  have  ensued. 
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CONCLUSIONS 

1.  Web-core  sandwich  type  cylindrical  pressure  hulls 
possess  structural  efficiencies  in  the  order  of  10-20% 
higher  than  those  of  conventional  ring-stiffened  construc- 
tion at  depths  of  2,000  feet  and  5%  higher  at  depths  of 
10,000  feet. 

2.  Higher  strengths  realized  with  the  sandwich  designs 
over  those  of  conventional  ring-stiffened  cylinders  can  be 
explained: 

a.  Symmetry  of  the  sandwich  cross  section  allows  a  more 
uniform  and  more  complete  stressing  of  the  available 
material.   The  outstanding  flange  of  a  T-frame  on  a 
conventional  ring-stiffened  pressure  hull  is  stressed 
only  in  the  circumferential  direction,  i.e.,  a  uniaxial 
stare  of  stress  exists.   In  a  sandwich  pressure  hull, 
the  "flange"  material  is  in  reality  a  second  shell 
which  is  more  efficiently  stressed  in  both  the  axial 
and  circumferential  directions,  i.e.,  a  biaxial  state 
of  stress  exists. 

2.  Sandwich  geometry  is  inherently  more  stable  structur- 
ally permitting  straining  of  the  material  well  into 
the  work-hardening  range  with  concomitant  beneficial 
effects  of  higher  strength  levels  (inelastic  materials). 
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3.  Inelastic  analysis  is  very  necessary  in  the  use  of 
high  strength  materials. 

4.  Sandwich  construction  permits  the  use  of  thinner  plat- 
ing material  with  its  superior  homogeneous,  higher  strength 
greater  ductility,  and  better  notch  toughness  character- 
istics over  thicker  plating.   In  addition,  thinner  plating 
is  easier  to  form  and  weld. 

5.  The  use  of  this  type  of  structure  is  appealing  in 
several  designs : 

(a)  Large  diamters  hulls  --  thinner  plating. 

(b)  Deep  depth  hulls  --  thinner  plating. 

(c)  Shallow  depth  hulls  (<2000  feet)  —  10-20%  more 
efficient  structure. 
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RECOMMENDATIONS 


1.  Utilize  the  membrane  stresses  developed  in  Appendix  A 
in  the  stress  analysis  and  failure  mode  programs  in  place 
of  the  outer  fiber  stresses. 

2.  Compare  experimental  results  to  theoretical  and 
determine  whether  membrane  stresses  or  outer  fiber  stresses 
are  more  realistic. 

3.  Formulate  a  more  sophisticated  optimization  program 
utilizing  a  search  technique  (reduce  program  cost) . 
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APPENDIX  A 

AXISYMMETRIC  STRESSES 
IN  AN  ISOTROPIC,  WEB  ■*- STIFFENED  SANDWICH  CYLINDER 
LOADED  WITH  UNIFORM  EXTERNAL  PRESSURE 


BASIC  EQUATIONS 
For  cylindrical  shells  loaded  under  uniform  pressure, 
the  three  principal  axes  of  strain  are 


x 


ax  "  v  (%  +  ar} 


<p         E  1  cp 


1 

E 


ar  '  v  (ax  +  Vj 


(1) 


Since  the  radial  stress  in  the  shell  is  generally 
small  in  comparison  to  the  longitudinal  and  circumferential 
stresses,  Equation  (1)  reduces  to 


x 


I 


va 


u =  M%  ■  voJ 


(2) 


for  the  longitudinal  and  circumferential  strains,  respect- 
ively. 
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The  equilibrium  of  forces  and  moment  of  forces  in  the 
element  requires  that 


v 


/  r 


T" 


dQ 


v 


N 


=  -  P 


dx    R 


dM 


x 


dx 


•x 


(3) 


if  the  beam-column  effect  due  to  the  axial  portion  of  the 
hydrostatic  pressure  is  neglected. 

The  moment  curvature  relationship 

M   -  -D  £fw 
x  dx/ 


where 


Et 


12JT=^Y 


(4) 
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Equations  (3)  and  (4)  can  be  combined  to  form  the 
following  general  differential  equation  of  equilibrium 

D^+^=P  ■      (5) 

dx^    K 

The  strains  in  equation  (2)  can  be  expressed  in  terms 

of  the  forces  N   and  N   as 

x       <p 

S  ■  re  %  '  VV  (6) 


The  circumferential  strain  can  also  be  expressed  as 


By  substituting  Equations  (6)  and  (7)  into  Equation 
(5) ,  the  general  equation  of  equilibrium  is 

D  +-—  w  =  P~5-Nv  (8) 

dx"    R2  R   x 

The  edge  conditions  at  X  =  0  and  X  =  L  are 

dw 1      ,  dwj       0 
dx|x=0    dx]x=L 


and 


-d3w        d3w I       Qx 
dx3 |x=0    dx3 |x=L   — 


With  these  boundary  conditions,  the  general  solution 
of  equation  (8)  is 
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w 


2R2Q 

wx 

EtL 

2 

(B1B3  +  B4  ■ 

"  B5  • 

"  B2B6 

r2p 

V 

N  R 

X 

Et 


E   t 


(10) 


where 


3(l-v2) 

2D2 


tzR 


B-,  = 


sinh  6  +  sin 


1    cosh  6  -  cos 


3^  = 


sinh  6  -  sin 
cosh  6  -  cos 


B-j  =  cos  (6— j  cosn  (Sr-J 


B4  =  sin  (eg  cosh  (eg 
B,  =  cos  (e=-J  sinh  (6r-) 


3.  =  sin  (e£)  sinh  (e£) 


STRESSES  AND  STRAINS 


The  longitudinal  stresses  on  the  shell  section  with 
externally  applied  pressure  for  the  outer  and  inner  surfaces 
are 
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From  Equations  (4)  and  (10) ,  tne  longitudinal  stresses 
of  Equation  (11)  are 


R  _x  -  ^x  1 
Cx    t  •  R     Up 

I 


> 


(12) 


From  Equations  (2),  (7),  (10),  and  (12),  the  circum- 
ferential stresses  are 


wnere 


R        Qx  (  Gl 

%  ^  P  +  L"  1G2  t   *    - 


*  L 


Gl    ■    6     (B233    "   B4    "   B5    +   31B6^ 
G2    =    6     (3^3    +    B4    -    B5    -    32B6) 


,=  /HZ^» 


(13) 


(14) 


DEFORMATION  OF  ISOTROPIC  CIRCULAR  RINGS 


For  a  web-stiffened  sandwich  cylinder,  the  outside 
and  inside  shells  are  separated  by  a  circular  ring.   A 
segment  for  this  type  of  ring  is  shown  below. 
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Y 
Y 


ring  segmcnt 


+ 


SMMl    EL£lA6NT 


Neglecting  stresses  in  the  X-direction  for  the  ring, 
Equation  (1)  reduces  to 


e  ,  =  — 


o  .  -   va 


£    =  ^7 


a   -  va 


^J 


(15) 


These  strains  can  also  be  expressed  as 


R 


£   = 


dw 
dr 


(16) 


The  equilibrium  of  radial  forces  on  the  small  element 
requires  that 

da 


a,  -  a  -   R-3 — 
cj)     r     dr 


=  0 


(17) 


With  the  relationships  in  Equations  (15)  and  (16) , 
the  general  solution  of  Equation  (17)  is 


W(R)  = 


T     fRRl]2     -i 
L"Qi  IRROJ  +  QoJ 


(1  -v) 


R 


RRO 


L 


-o    fRRii  +  n    fRRIT 

Qi  [RROJ     +  Q0  [rROJ 


(l-v) 


rRRI 
I  R 


2   RRO 


1  - 


fRRll2"1 

[rroj 


(18) 


where 


Z  = 


,RRO 

RRI 

RR 


1-v2 
RO  + 

RI  - 

RRI 
RRO 


TO 
2 

TI 


By  differentiating  Equation  (18)  and  by  combining 
Equations  (15) ,  (16) ,  and  (17) ,  the  circumferential  stress 


is 


■Q±(RR)Z+QQ 


VR) 


R  ! 


r 


RRO 


J  T  L"Qi 


Q, (RR)+Q  (RR) 
o 


Trri 

R 


bf  R  1 
2 [RROj 


1-(RR): 


(19) 


and  the  radial  stress  is 
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ar(R) 


-Qi(RR)'-:-Q0 

L 


[rroJ 


■Qi(RR)+QQ(RR) 


R 


RRO 


l-(RR) 


RRI 


R 


(20) 


DEFORMATIONS  OF  SANDWICH  SHELLS 

3efore  determining  the  stresses  in  the  cylindrical 
shells  of  a  sanawicn  structure,  it  is  first  necessary  to  ob- 
tain the  forces  Q   and  X   on  the  outside  shell  and  the 

o      o 


forces  Q.  and  N.  on  the  inside  shell. 


OuUi^ie    S^ell 


-o.   x 


jj     minn  _l/_ 


W«k> 


"~1 


L 


f 


ni  r 


_,    —   j—  — 
-  i. 

r 


_Ai 


I 


L 


-„* 


"R. 


"R: 


The  radial  deflection  of  the  outside  shell  at  the  web- 
shell  juncture  can  be  expressed  interms  of  the  forces  as 


W   =  g  Q   +dN   +  f  P 
os    ^o  o    o  o    o 


(21) 
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and  the  radial  deflection  of  the  inside  shell 


W. 

is 


g-Q- 
3i  i 


d.N. 
1  i 


(22) 


Similarly,  the  radial  deflection  of  the  outer  fiber 
of  the  web  at  the  web-shell  juncture  can  be  expressed  as 


W 


ow 


=  g 


oo 


L 


o 


+  fe  P 


+  g  -Q • 


(23) 


and  the  radial  deflection  of  the  inner  fiber  can  be 
expressed  as 

W .  -a.         Q   +  5-  P   +  g .  .  Q . 


(24) 


Since  the  stresses  in  a  sandwich  structure  are  linear 
functions  of  the  applied  pressure,  the  coefficients  of  the 
force  in  Equations  (21) ,  (22) ,     (23),  and  (24)  will  be 
developed  for  a  unit  ©f  applied  pressure  (p  =  1.0) .   By 
substituting  the  geometric  properties  of  the  outside  shell 
into  Equation  (10) ,  the  coefficients  of  the  forces  in 
Equation  (21)  are 


(25: 


So  = 

-  2Rcreo 

sinhe  -i-  sinG 
o        o 

2ETOL 

cosh 8  -  cos  6 
o       o_ 

d  = 
o 

-  vRO 
ETO 

f   = 

RO2 

o 

ETO 

where 


=  it 


o    /3(l-v^)  /lO^TO 
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Also  by  the  geometric  properties  of  the  inside  shell, 
the  coefficients  of  the  forces  in  Equation  (22)  are 


-  2RI26. 


*i  = 


2ETIL 

vRI 

ETI 


sinh0.  +  sin9. 
l l 

coshB.  -  cosG. 
1.       i 


where 


6i  =  /3(l-v2)   /Ri.TI 

The  hyberbolic  and  trigonometric  function  in 
Equations  (25)  and  (26)  are  redefined  by 


F  = 


sinhe  +  sini 


coshe  -  cosi 


(26) 


(27) 


By  Equation  (27),  the  coefficients  g   and  g.  are 


*i  = 


-2R02 
ETOL   O 


■ml  F 

ETIL   i 


(28) 


By  substituting  combinations  of  unit  and  zero  forces 
into  Equation  (17)  ,  the  coefficients  in  Equations  (23)  and 

(24)  are 

g    =  CRRO  [(1  +  RR2)  -  v(l  -  RR2)] 
oo 


gQi  =  -2CRRI(RR) 

gio  =  2CRR0(RR) 

g±   =  -CRRI  [(1  +  RR2)  -  v(l  -  RR2)] 


(29) 


§1 

where 

U    bE(l  "  RR2) 

From  the  equilibrium  conditions  of  the  end  forces 
on  a  sandwich  cylinder  loaded  under  an  external  uniform 
pressure,  the  relationship  between  the  outside  and  inside 
axial  forces  for  p  -   1.0  psi  is 

on  2 

N  R   +  N.R.  =  ^-  (30) 

o  o     i  i     2 

By  assuming  that  the  axial  membrane  strain  in  the 

outside  shell  is  equal  to  the  axial  membrane  strain  in 

the  inside  shell,  i.e.,  r~K  -   ^t^f/    the  axial  forces  are 

ETO    ETI 


N_  = 


RO 
2 


O         TI aRI 

to'ro 


TI 
N.  =   —  N 
i    TO  o 


(31) 


For  compatibility  of  deformations  in  the  radial  direc- 
tion, the  relationships  that  must  exist  are  W   =  W  ,  and 

os     ow 

W.   =  W.  .   With  these  relationsnips  and  Equations  (21)  to 
is     iw  ^ 

(29),  the  radial  forces  Q   and  Q.  can  be  expressed  in  terms 
of  the  two  following  simultaneous  equations: 


86. 


(g  -g  )  Q  -  g  -Q-  = 


-RO 


ETO 


1  - 


v/2 


1  + 


TI.RI 

to'ro 


+  77  g 
2  yoo 


v  RO' 


-9ioQo  +  ^i^ii^i  =  2E  fl 


1  - 


1  + 


TI  RI 

to'ro 


b 

2  ^io 


+  -*  g 


(32) 


The  forces  N  ,  N . ,    Q  ,  and  Q.  can  be  computed  from 
Equations  (31)  and  (32) . 


STRESS  DISTRIBUTION  IN  SANDWICH  STRUCTURES 


The  stresses  on  the  outside  shell  of  the  isotropic 
sandwich  are  determined  by  substituting  N  ,  Q  ,  and  the 
geometric  properties  of  the  outside  shell  into  Equations 
(12)  and  (13).   From  these  relationships,  the  circumfer- 
ential stress  in  the  outside  shell  on  an  outer  fiber  is 
a 


4oo  -  -  P  i  £  -1  +  IT  <S20  +  f  G10>  1 


on  an  inner  fiber  is 


RO 


%oi  =  "  ?  *§  [  -1  +  L 


and  the  membrane  stress  is 

RO 


o 


(G 


20 


f  G10}  3 


(33) 


<J)om 


p  ¥§  I  -1  +  TT  G20  > 
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The  longitudinal  stress  in  the  outside  shell  on  an 
outer  fiber  is 


_  RO    N     Q   _ 

axoo  =  P  TO    -£-  ^°-G10 

[_  RO    n^  J 


on  an  inner  fiber 


xoi 


=  P 


RO 
TO 


N     Q 

-2.  +  -2.  r 
RO    Lip  ^10 


and  the  membrane  stress  is 


(34) 


N 

a    =  P  -2. 
xom       TO 


And  likewise  the  stresses  on  the  inside  shell  are 
determined  by  substituting  N.,  Q. ,    and  the  geometric 
properties  of  the  inside  shell  into  Equations  (12)  and  (13) 


P    St  ?i 

°$io  TI    L 


G„.    +    7-  G,  . 

2i         \p      li  _ 


<|>ii 


Qi  r 

-  P   £=•  ~    I    GnJ    -  £  G-,  , 


RI    "I 
TI    L 


ip   "li 


c})im 


-  p  £L?i  g 

r    TI     L     ^2i 


(35) 


xio 


Xll 


a 
xim 


■   P 


=   P 


=   P 


RI 


Ni        Qi 


TI      _RI         Lip      li 


RI 
TI 


N. 

l 

TI 


Ni 

RI   +   Lip   Gli 


(36) 
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The  stresses  in  the  web  stiffeners  at  the  outside 

and  inside  shell  locations  are  determined  by  substituting 

the  outside  radius  RRO  and  inside  radius  RRI  into 

Equations  (19)  and  (20) .   Also,  the  total  force  on  the 

outer  surface  of  the  web  consists  of  the  radial  force  Q 

o 

at  the  web-shell  juncture  and  the  force  attributed  to  the 

applied  pressure  directly  on  the  web.   The  total  forces 

(H   and  H.)  on  the  outer  and  inner  surfaces  of  the  web 
o      i 

for  a  unit  pressure  (p  =  1.0  psi)  are 

H   =  Q   +  b/2 
o   vo    ' 

Hi  =  Qi  (37) 

With  these  relationships,  the  circumferential  stress 
in  the  web  on  the  outer  surface  is 


$wo 


(-Hi(RR) 2+HQj  +  (-Hi(RR)+HQ(RR)) (RR) 


|  (1-  RR2) 


and  the  radial  stress  (33) 

(-Hi(RR)2  +HQ)  -  (-Hi(RR)+HQ(RR))  (RR) 


a_  =  -P 


rwo  b 

2  Cl-RR2) 


and  on   the   inner   surface 


<J>wi 


(-Hi(RR)2+HQ)  (RR)     -     (-Hi(RR)+Ho(RR)) 


j    (RR)      (1-RR2) 
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rwi 


(-Hi (RR) 2 +Hq) (RR)     -     (-Hi(RR)+HQ(RR)) 


b 


(RR)      (1-RR2) 


(39) 
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APPENDIX  B 


BUCKLING  ANALYSIS 


1.  inelastic  Asymmetric  (Lobar)  Buckling  Analysis. 

The  buckling  equations  for  a  fully  plastic  cylinder 
specialized  for  the  case  of  hydrostatic  pressure  loading 


^T    [l    8rJ_u         1    32v 
2E0  I  2    3x2         r   3x3<}) 


32u 


1_  _3_w 
R    3x 


1     32v 


,    3    3'u  1 

4   3x2       JbJ  jpr  T   4R  5x9(|j 


=    0 


(la) 


E 


T    (l    32u 


1     3zv 


2,-, 


1     3w 


E       [2    3x3<j>         R2    3£ 


2    +    R2    "3<j>J 


1    c^v  1     32u 


4    3x2         4R    3x3<J> 


(lb) 


3      R 


fl    3u        1    3v        w 
2    3x        R    34)         R 


+   D 


E       [4    3X1*         R23-x234)2    +   r*    a^"( 


+  ii^- 


4    3X1*         R2    3x23cj)2 


+   N 


J 


32w 
x    3x2 


dc) 


+    %    57  f^   +    P    =    ° 
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where  x  and  4>  are  respectively  -he  axial  and  circum- 
ferential coordinates, 

u,  v,  and  w  are  the  axial,  tangential,  and  radial 
di  splacements , 

E   and  ET  are  the  secant  and  tangent  moduli, 

R  is  the  radius  to  shell  mid-surface, 

t  is  the  shell  thickness, 

v  is  Poisson's  ratio, 

D  is  the  bending  rigidity  =  E  t3/12(l-v2) 

9 


N   and  N^  are  forces  per  unit  length  in  the  axial  and 
circumferential  directions,  and 

p  is  the  hydrostatic  pressure. 


With  several  differentiations,  u  and  v  can  be  elimin- 
ated from  Equations  (la) ,  (lb)  and  (lc)  so  that  a  single 
eight-order  equation  in  w  is  obtained: 


D 


E  f 

=£  V8w   +      1- 


E. 


^f3   3*w       J^_     a^w 
1 V     [2    dxF        R2    3x23cj>2_ 


*i 


E 


-1 


'T 


(3    3 


'w 


1 


38w 


4    3x8    T  R2    3x6 3(p 


E    t 
s 

R2" 


3"w 


+   N 


+    N 


V      3x2    +    4 


rE 


E, 


-1 


1 


32w 


R^     3cp2  4    R2 


[ET 


36w 
3x6 


36w 


3x43cp 


=    0 


(2) 
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where  V       indicates 


1      8 


3x2         R2    3<J>2 


A  solution   to    this    equation   can  be  written: 

W  =   A   sin   K<j>   sin    Xx  (3) 

where        K  =  N 

,         mir 
A   -   — 

L  is  the  length  of  the  shell,  and 

m  and  n  are  integers . 

This  solution  satisfies  the  conditions  of  simple  sup- 

32w 


port  at  the  ends  of  the  cylinder,  i.e.,  that  w  and 
vanish  at  x  =  0  and  x  =  1. 


3x' 


These  conditions  are  not  unreasonable  for  stiffened 
cylinders  since  it  is  likely  that  the  effective  rotational 
restraint  will  be  limited  by  the  formation  of  plastic 
regions  arising  from  high  bending  stresses  near  stiffeners 

By  substituting  the  solution  (3)  back  into  (2) ,  the 
following  characteristic-value  equation  is  obtained: 


+  a2!  + 


1- 


'k2 

R2 


+  X 


f3X2    K2)   3X1* 


L£  -1 
[*T    J 


3X 


K 


2  "I 


4   +  R2 
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E  t 
+  _s_  X*  -  N 
R2 


R2 


1  2 

+  A2   A2  + 


E 


-  1 


'T 


N   R2 


To  simplify: 


R* 


+  A 


+  |  -S.  -  l' 

4  IE 


=  0 


(4) 


$  E 


A2  + 


K2 
R7 


1  + 


m2lT2R2 


f  E 


c  = 


fEs 
Em 


-  1 


(5) 


Now  since: 


Nx  "  'x'  =  KxPt 


the  equation  is  then  rearranged  so  that  an  expression  for 
P  ,  the  plastic  buckling  pressure,  is  obtained: 


2    ET 
P      =    fDA2    ~ 

P  Eo 


1  +  c#   11  +  *  +  4 


[ 


1  - 


E    tfS* 
s 

R2A2 


Kxt* 


1    -    $     (l    - 


.1  r 

f)JL1  + 


C$' 


(6) 


Equation  (6)  can  now  be  minimized  for  P  with  respect 
to  N/M  by  setting 
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**     = 


m***   ET    f/Rtl1* 

^P-  -   o 

"  1   -    2*       (1-f )    +   ^c$    2 
P                        4      p 

9        2S     t     LJ 

3-2$    (1-f)    +   ic$    2 
p                     4      p 

(7) 


(8) 


and 


P      = 
P 


4m27T2ETf 
9  3> 


'  t  ' 

Rt 

RK 
I       x, 

L 

C  9  o 

1   +   -J£    (3   +    $      +   |c$    2) 
4        p         4      p 

3-2$       (1-f)    +   |c$    2 
p  4      p 


(9) 


Since  P   in  (9)  is  proportional  to  m2 ,   m  must  be 
equal  to  one  in  all  cases  when  N  is  greater  than  0. 

After  some  investigation,  an  approximate  value  of  <p 
was  determined 


9    -    1.23 


/Rt 


(10) 


Equation  (9)  can  be  further  simplified  to 


P  = 
P 


'"2V  [  t  ] 
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RK 


x 


/Rt 


2  r~ 


1  + 


3c$   -, 


3  -  2$(l-f) 


(11) 
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A  similar  analysis  for  the  elastic  region  yields: 

7T2Ef  f    t    )      f/Rtl2     f  1 

1-V*)«    l^U     I   L    J        l3-2*(l-f)J 


?e  ' 


(12) 


Although  Equations  (11) ,  (12)  define  the  buckling 
pressure  for  the  plastic  and  elastic  regions,  no  solution 
is  given  for  the  inelastic  region  which  lies  between  these 
two  limiting  cases.  However,  by  employing  an  empirical 
correction  factor  wherein  Poisson's  ratio  is  regarded  as  a 
variable,  one  can  arrive  at  an  expression  which  reduces  to 
the  proper  limiting  values. 


v  = 


2     ej 


(13) 


Thus  rewriting  Equation  (11) , 


P      = 


7T^ETf 

■3(l-v2)     [ 


RK. 


x 


f/RtV 


3£ 
1+4 


r5^-i 


3    -    2*(l-£) 


(14) 


2 •  Inelastic  Axisymmetric  Buckling  Analysis' 
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The  differential  equations  of  equilibrium  describing 
the  plastic  buckling  of  a  cylindrical  shell  subjected  to 
external  hydrostatic  pressure  are  expressed  as  Equations 
(1). 

Equations  (1)  admits  to  a  solution  of  the  form 


u  =  A   cos   N4>  cos    Xx 

v  =   B   sin  N#  sin    Xx 

w  =   c   cos   N<j)  sin   Xx 
where    X   =   Mtt/L,      n   an   integer 


(15) 


This  solution  describes  a  buckling  mode  with  M 
halfwaves  along  the  cylinder  axis  and  2N  halfwaves 
around  its  circumference.   It  also  satisfies  boundary 
conditions  of  simple  support  at  the  ends  of  the  cylinder; 
i.e.,  w  and  a2 w/ ox2  vanish  at  x  =  0  and  x  =  L. 

The  buckling  condition  of  the  shell  is  represented 
in  Equation  (4) . 

For  axisymmetric  buckling  N  =  0  and  for  m  =  1,  the 

minimum  value   of    X   is    X  =   tt/L. 

mm    ' 

Substituting  this  value  into  Equation  (4)  results  in 


ETtL2 
N      = 


9tt2R2 


2.    2„h  i  3 


R't'tt»(l    +    go)  jL{1+6) 


L*  (1  +   |c) 


(16) 
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replacing   N     with   a    t  =   K  pt 


p     = 


ETL' 


9K    7T2R2 

X 


'VR't'g  +  jo      9(1  +  c)- 


Again  using  the  expression 


(1  +  Jc) 


(17) 


1    E_ 
V  =  2  "  E"   2 


s   1 


EmL' 


Po  = 


12K  (1-V2)7I2R2 


ir'R'tM-I-Hfc) 


12(l-v2) (1+c) 


(18) 


3.  Web  Buckling  Under  Eydrostatic  Loading* 

Two  modes  of  buckling  will  be  considered: 
(a)  Buckling  of  a  Ring 


N-    Z 
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PC  "  (n  _1)   T^T  " 


LR; 


m 


(19) 


(b)  Web  Buckling  under  Edge  Compression 


cry 


%    "^ 


LU  11  Hill         _i. 


t  t  t  fTTTTTT 

i a  » 


; 


0 


I- 


a  =  z<fR 


w 


o     m2    +  a      n2    hr 
x  y  Id 


2  „2 


a2b 


2 

d7 


(20) 


where 


Therefore 


o      =0,m>l,n=l 
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a      = 
x 


7T2a2D    (m2        n 


m2b        a' 


2 

dr 


4tt2D 
d2b 


(21) 


a      =0,111=1,11  =   1 


=  lzJlR  f«ii  +  ail  -  zi£ 

d2 


n*b        a' 


d2b 


(22) 


X  y  d2b 


(23) 


and 


thus 


a      =   cj  .    =  K. 
x  cp  cj)p 


a      =   a     =  K  _ 
y  r  rP 


P(K,    +    4K    )    = 


4tt2D 
d2b 


P      = 


4tt2D 
d~?b~ 

(Ka    +    4Kr) 


(24) 


where 

d  =  h 

Eb3 


D  = 


thus 


12(l-v2) 


*Reference  (2),  (4),  (8),  and  (22). 
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P      =  4TT2E(b2/h2)  (25) 

C         12     (1    -    v2)     (K      +    4Kr) 
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APPENDIX  C 

MATERIAL  PROPERTIES* 

The  properties  for  the  four  materials  utilized  in  the 
optimization  studies  are  tabled  below.   The  secant  and  tan- 
gent moduli  are  defined  as 

E  -!i 

s    e . 

i 


do 


Em  = 


T    de 


*i 
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(1)  HIGH  STRENGTH  STEEL  (a   =  160,000  psi) 
a .  /a 

1    y 

0.65 
0.70 
0.75 
0.80       1.00       1.00        E  =  30  x  10*  psi 


ET/E 

/EsET/E2 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.95 

0.98 

0.76 

0.86 

0.55 

0.72 

0.30 

0.52 

0.12 

0.30 

0.06 

0.18 

0.85       0.95       0.98        \>  =    0.30 

0.90       0.76       0.86        p  =  490  pcf 

0.95 

1.00 

1.05 

1.10 


(2)  TITANIUM  ALLOY  (    =  110,000  psi) 


1  ■    i 
ET/E 

1.00 

a  •  /a 

1  y 

/EsET/E2. 

0.65 

1.00 

0.70 

1.00 

1.00 

0.75 

0.85 

0.98 

0.80 

0.68 

0.85 

0.85 

0.55 

0.75 

0.90 

0.44 

0.62 

0.95 

0.32 

0.50 

1.00 

0.22 

0.41 

1.05 

0.15 

0.30 

1.10 

0.06 

0.22 

E  =  18  x  105  psi 


v„  =  0.30 

hi 


276  pcf 
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(3)  ALUMINUM  ALLOY  (a   =  65,000) 

CTi/gy 

0.65 
0.70 
0.75 
0.80  1.00        E  =  10.8  x  iob  psi 


ET/E 

/EsET/E2 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.62 

0.80 

0.34 

0.55 

0.14 

0.28 

0.05 

0.12 

0.03 

0.06 

0.85  1.00       v„  =  0.30 

E 

0.90        u*°*       0.80        p  =  173  pcf 

0.95 

1.00 

1.05 

1.10 


(4)  HY80  STEEL  (a,  =  80,000  psi) 


ET/E 

1.00 

ai/Gy 

/EsET/E2 

0.65 

1.00 

1.00 

1.00 

1.00 

1.05 

0.00 

0.00 

E  =  30  x  106  psi 


v„  =  0.30 


490  pcf 
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APPENDIX  D 


SAMPLE  OUTPUT  OF  COMPUTER  PROGRAMS 


1.  STAS  -  Stress  Analysis  Program 

2.  FMCS  -  Failure  Mode  Criteria  Program 

3.  OPTI  -  Optimization  Program 
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STRESS    ANALYSIS    OF    A   WEB    STIFFENED    SANDWICH  CYLINDRICAL    PRESSURE    HL 

OUTSIDE  RADIUS* IN. 1=      61.50 

INSIDE  RADlUSUN>«      58.50 

OUTSIDE  THICKNESS  (1*0*       0.60 

INSIDE  THICKNESS*  0.60 

WEB  THICKNESSCIN)*       0.60 

WEB  SPACIN&(IN)>  5.40 

MXOO    MPHIOO    MXOI    MPHIOI    NxOl    MPHIIO    MX  1 1    NPHUl 

-37.78   -47.87   -14.75   -40.96   -16.95   -4i.9i   -35«,S8   -47.30 
FXOO    FPHIOO    Fxpi    FPHIOI    FXlO   FPHIIO    Fxll    FpHIIl 

-3.19   -36.45   -49.35   -50.30   -44.95   -51.19    -7.58   -39,98 
HRW       HS       HTS       BTS       IEFF       D 

0.050  0.500  5.000  1.000  16.9  3*600 

WB«  0.377 

KR=  3.071  KP=  83.887 

STRESS  ANALYSIS  OF  A  WEB  STIFFENEO  SANDWICH  CYLINDRICAL  PRESSURE  HL 

OUTSIDE  RADIUS! IN. )»  123.00 

INSIDE  RADIUSUN)*  117.00 

OUTSIDE  THICKNESS  IIN1  =       1*20 

INSIDE  THICKNESS-  1.20 

WER  THICKNESSIINI*       1.20 

WEB  SPACINGUN)*  10.80 

MXOO    MPHIOO    HXOI    MPHIOI    MXOl    MPHIIO    MXII    MPHIII 

-37.78   -47.87   -14.75   -40.96   -16.95   -41.91   -35.58   -47.50 
FXOO    FPHIOO    FXOI    FPHIOI    FxlO   FPHI 10    FXII    FPHIH 

-3.19   -36.45   -49.35   -50.30   -44.95   -5U19    -7.58   -39.98 
HRW       HS       HTS       BTS       IEFF       D 

0.050      0.500      5.000      1.000      270.6      7.200 
WB=      0.377 

KR«      3.071       KP«     83.887 

STRESS  ANALYSIS  OF  A  WEB  STIFFENED  SANOWICH  CYLINDRICAL  PRESSURE  HL 

OUTSIDE  RADIUS! IN.»»     184.50 

INSIDE  RADIUSlIN)-     175.50 

OUTSIDE  THICKNESS  (IN)*       1.80 
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INSIDE    THICKNESS*  1.80 

WEB  THICKNESS(IN)*  1.80 

WEB  SPACINGUNI*      16.20 

MXOO    MPHIOO    MXOI    MPHIOI    MXOI    MPH1I0    NXII    MPHI 1 1 

-37.78   -47.87   -14.75   -40.96   -16.95   -41.91   -35*58   -47.50 
FXOO    FPHIOO    EXOI    FPHIOI    FxIO   FPHIIO    FXII    FPHIII 

-3.19   -36.45   -49.35   -50.30   -44.95   -51.19    -7.58   -39.98 
HRH       HS       HTS       8TS       IEFF       0 

0.050      0.500      5.000      1.000     1369.9     10*800 

W3=      0.377 

KR*  3.071  KP=  83.887 

COMPILE    TIME*  7.53    SEC, EXECUT ION    TIME*  0.86    SEC. OBJECT    C00E«  9856    BVTI 
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